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SUMMARY 


Equations of motion are derived for free vibration of a uniform straight 
free-free beam with a rotational mass attached at each end. The mode-shapes are 
shown to be linear combinations of trigonometric and hyperbolic sine and cosine 
functions. By using boundary conditions, nonlinear algebraic equations are derived 
for both lateral and torsional vibrations. These solutions can then be used to 
compute natural frequencies and mode shapes. 

A computer program was written to solve the nonlinear algebraic equations. 
This program uses the Newton-Raphson method to compute natural frequencies, mode 
shapes and node points for any number of modes, for any given parameters. 


INTRODUCTION 

The lateral and torsional vibration characteristics in the design of certain 
systems subject to external Interference are not fully understood. Some common 
examples of such systems are the wing of an aircraft, the base of reciprocating 
engine, and the hull of a ship. Some of these systems can be effectively 
approximated by a single beam with varying end conditions. This paper will deal 
with the characteristics of free vibration of such beams. 

Both the properties of natural frequency and location of node points 
provide the designer with those tools needed to either minimize or maximize the 
vibration of a system. Knowledge of the natural frequency of a system allows a 
designer to avoid resonance. If a concentrated force is located at a node point, 
then the induced vibration will be minimized. These characteristics are easily 
obtainable from tables of predetermined eigenvalues for simple beams such as 
cantilever, hinged-hinged, and a few others (refs. 1, 2, 3). However, the problem 
is more complicated for those configurations having something other than a hinged, 
clamped, or free end mass. 

This paper describes the derivation of the differential equations for free 
vibration. The linearized equations satisfying the various boundary conditions, 
will then be obtained using the methods discussed in reference 1. The 
configuration to be considered is a free-free beam with rotating masses at both 
extremities. This configuration is considered to be the most general case, with 

all the classic cases able to be represented by the lower and upper limiting 
cases. The computer program FFBEAM which computes mode shapes, nodal points and 
natural frequencies of free vibration is discussed in the text of this paper. A 
listing of FFBEAM and sample input and output are also included in the appendices. 
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cross sectional area of beam or shaft 
coefficients of mode shape equation 

dimensionless position variable 
stiffness of beam 
modulus of rigidity 
polar moment of inertia 

inertia of mass at x = 0 

inertia of mass at x = L 

length of beam or shaft 
mass at x = 0 on beam 

mass at x = L on beam 

moment 

moment of x = 0 on beam 

moment of x = L on beam 

time solution of partial D.E. 

lateral mode shape 

time 

torque 

torque at x = 0 on beam 

torque at x = L on beam 

lateral displacement of beam 
shear force 

shear at x = 0 on beam 

shear at x = L on beam 

location on axis of beam or shaft 
eigenvalue matrix 
phase angle (rad) 
eigenvalue variable 
dimensionless position variable 
torsional mode shape 
density of beam or shaft 
angular displacement of shaft 
rotational velocity (rad/sec) 
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DIFFERENTIAL EQUATION OF LATERAL FREE VIBRATION 


The basic differential equation of lateral free vibration is best discussed 
using figure 1. Therein, a section of beam of length dx can be seen. Its 
longitudinal displacements and slopes are assumed to be negligible. Then, setting 
the transverse shear force acting on the element equal to the element's mass times 
acceleration, the following partial differential equation is obtained (ref. 1). 


9 v(x,t) 

3t2 


PA 


5 v(x,t) 


( 1 ) 


where v(x,t) is the lateral displacement. 

This equation is solved by separation of variables. That is, v(x,t) can be 
written as (ref. 1): 


v(x,t) = p(t) r(x) 


( 2 ) 


Then substituting the right half of equation (2) into the original D.E. and 
rearranging the variables. 


d^p(t) 1 _ ^ 1 d^r(x) 

dt2 ^ PA r(x) ^^4 


But (3) can be true if, and only if, each side is Identically equal to a constant. 
If this constant is -o)2^ ordinary homogeneous equations are obtained: 


and 


= 0 
dt 


- 6^r(K) = 0 
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The solutions for (4) and (5) are (ref. 1): 


r(x) = AjSin 3x + B^cos Bx + C^sinh Bx + D^cosh Bx 


and 


p(t) = cos(o)t - a) 


( 7 ) 


where a represents the phase angle. But this solution is more convenient to use 
in its nondimensionalized form. Therefore, the new variable e = x/L, where L is 
the length of the beam, is introduced. The solution of displacement after the 
substitution is found to be 


r(e) = AjSin + Bj^ cos + C^sinh B^^ + D^^cosh 


( 8 ) 


where 


B 


4 

2 


pA 

El 


2 

(jj L 


(9) 


Boundary Conditions for Lateral Vibration 


The configuration being considered is a beam with rotational masses at both 
ends. A sketch of the configuration is shown in figure 2. There the two free body 
diagrams show the anticipated reactions at the masses. 

Since the original differential equation is of fourth order, four boundary 
conditions are needed. By using the equations for shear and bending moment 

(ref. 1): 


( 10 ) 


( 11 ) 


one can obtain the four necessary linear equations. 


El d v(£,t) 




M = 


El d v(c,t) 


de' 
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I 


Following reference 1 and referring to figure 2, the first boundary condition 
involves the shear Vi: 


„ _ El d^v(e,t) d^v(e,t) 

V, - - — ;; - m, r at 


( 12 ) 


de' 


dt 


By substituting the right-hand side of equation (2) into equation (12), the 
following is obtained: 


,3 


d^r(e) 

de^ 


p(t) = m^r(e) 


d^p(t) 

dt^ 


(13) 


Equation (4) can be rewritten as: 


d^p(t) 
dt^ - 

P(t) 


(14) 


Using equation (14) in equation (13), we have: 


2 / V El d r(e) 

03 m r(e) - ^ = 0 


or 




where &2 given by equation (9). 


(15) 


(16) 


II. The derivation of the second boundary conditions employs the same procedure as 
that of the first boundary condition. The only difference is that they are 
applied to the mass m£ and shear V 2 . The resulting equation is: 


d^r(c) 

de^ 


/ 

■d 


at 


(17) 
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III. Following reference 1, the third boundary condition involves the moment 
(see figure 2): 


= I^0(e,t) 


,1 


d^r(e) 

dc2 


p(t) 


at e = 0 


( 18 ) 


By assuming small angular displacements, one may represent 0(e,t) by the slope 

ft 

3v(x,t)/8x. Therefore, using equation (2), for small displacements, 6 can 
be approximated by: 


dr(x) 

dx 


p(t) 


thereby obtaining: 




dr(£) d^p(t) _ El d^r(e) 

de J.2 ,2 ,2 

dt L de 


p(t) 


(19) 


By using equation (14), equation (19) can be reduced: 


or 


d^r(e) 


de 


d^r(e) 


de‘ 


Ij L m 
El 


dr(e) 

de 



dr(e) 

de 


at e = 0 


( 20 ) 


( 21 ) 


IV. The fourth boundary condition uses the same method as the third. After 
employing the same steps as outlined above, the following is obtained: 


d^r(e) 


de" 



dr(e) 

de 


at e 


1 . 


( 22 ) 


By substituting equation (8) and appropriate values of e (i.e., e = 0 at one 
boundary and e = 1 at the other boundary) into equations (16), (17), (21) 

and (22), the following four linear equations are obtained: 
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( 23 ) 
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(25) 


1202 \ f h4 

^ CO. 6,- Sin 6,jA, .1^ sin 6,- co. 6,|B^ 


^2®2 ] I 4h 

+ 1 j cosh g^+ sinh 02/^1 "'' 1 3 ^inh cosh ^ 2)^1 ~ ® 


PAL' 


PAL' 


(26) 


Obtaining Nontrivial Solutions 


Equations (23)-(26) can be compactly written in the vector-matrix form as 
follows: 


Z(g 2 ) 


l“i 


= 0 


(27) 


where Z(g) is the 4x4 coefficient matrix whose entries are functions of g. 
Non-zero solutions (Aj , Bj , Cj , D^,) exist only when the determinant of 
Z(g) is zero. Therefore, the first step in obtaining non-trivial solutions is to 
obtain the real solutions of the nonlinear equation: 


det 


_z(g2) 


0 


(28) 


where det [ ] denotes the determinant. 
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A solution B* is substituted back in equation (27), and a degenerate system 
(usually of rank 3) of algebraic equations is obtained. Making Aj = 1, the 
remaining three coefficients B][ , and Dj can be uniquely determined 

for each B*. The computer program MACSYMA was used to carry out these proce- 
dures. A listing and description of the steps can be found in Appendix C. The 
resulting determinant of Z(B) obtained in MACSYMA is used in the program FFBEAM 
which is discussed below. 


PROGRAM DISCUSSION 

The computer program FFBEAM (see Appendix A for a listing) was written for the 
purpose of calculating non-trivial solutions using the method described in the 
preceding section. FFBEAM calculates the frequencies of vibration, mode shapes and 
the locations of nodes for any number of modes. 

The program uses a simple root- finding subprogram (NONLIN) that employs the 
Newton- Raphson method to find the eigenvalues of the determinant. The size of the 
equation representing the determinant necessitated the separation of the equation 
into three separate quantities that are subsequently added together. FFBEAM uses 
the subprogram GELIM (for Gauss (-Seidel) elimination) to solve the degenerate 
system of equations (eq. 28) discussed in the previous section. 

After computing the values of B, A^ , Bj , C^ and Dj from 
equations (27) and (28), the program uses simple substitution to find the 
frequencies. The subprogram NONLIN is again used to find the node points of the 
mode shape equation. A further explanation of the input as well as output is given 
below. 


PROGRAM OPERATION 

Input .- A total of eight input variables is needed to run FFBEAM. A 

consistent set of units is needed for the input (meters-kilograms-seconds or 
f eet-pounds-seconds is recommended) • The necessary input information includes two 
concentrated endpoint masses (Ml and M2), two endpoint inertias ( II and 12), 
beam mass and length (MU and L, respectively), beam stiffness (El), and the 

number of modes solutions are needed for (NMODE). Because of the small number of 
input variables the input can simply be written in the text of the program. A 
listing of sample input is given in Appendix B. 

Output .- The output computed by FFBEAM includes the natural frequency of 
vibration, eigenvalue, mode shape equation and the location of node points for each 

mode. A listing of sample output is given in Appendix B. The first three pages of 

Appendix B include the frequencies, node positions, eigenvalues and mode shape 
equations. The final pages of Appendix B are the plots corresponding to each of 
the modes. These plots were found useful in identifying the relative motion of 
points of interest on the beam. 
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THE DIFFERENTIAL EQUATION OF TORSIONAL VIBRATION 

The differential equation of torsional vibration will be derived next, using 
the method described in reference 1. This method utilizes a differential element 
of a shaft of length dx (as seen in fig. 4) and the dynamic equilibrium of this 
element. 


It is known that the torque at any position of a circular shaft may be written 
as (ref. 1); 


T = GI 


3<l>(x,t) 


9x 


(29) 


Utilizing equation (29) and the differential element shown in figure 3, it follows 
that the net torque acting on the element at any time is (ref. 1); 


|1 ax.oi 

9x 


(30) 


This torque is opposed by an inertial torque: 


- I p 


(31) 


Hence, sum of these torques must be zero, and the differential equation governing 
torsional vibration of circular shafts is: 


9^(|>(x,t) _ G 9^([)(x,t) 

9t2 


9x2 


(32) 


This equation can also be solved by separation of variables. The angular 
displacement of the shaft is written as: 


4>(x,t) = p(t) 9(x) 


(33) 


Then, after substituting into equation (32) and rearranging, gives: 

1 d^p( t) ^ 1 d2s(x) 

p(t) ,2 p 0(x) 2 (34) 

at dx 

The solutions for both 0(x) and p(t) are obtained by setting both sides 
equal to -oj2; 

p(t) = cos(wt - a) ( 35 j 

0(x) = Aj sin X + cos m ^ x (36) 
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Using the dimensionless variable e = x/L, equation (35) can be rewritten as: 


0(e) = sin 0^^ 6^^ 

where 



( 37 ) 


(38) 


Boundary Conditions of Torsional Vibrations 

The configuration being considered, as seen in figure A, is a shaft with a 
rotational mass at each end. The method used to obtain the two linear equations is 
outlined in reference 1. The only relationship needed to obtain the boundary 
conditions for this system is (ref. 1): 


GI 


T = 


d6(g) 

de 


= - I to"^6(e) 


(39) 


I. The first boundary condition involves the mass located at e = 0 (that is, at 
the left end). Using the relationship in equation (39) one obtains: 

GI 

d9^ ^ ^ q 

1 L uG 1 

which can be rewritten as 


de(£) 

de 




I L p 
P 


at e = 0 


(41) 


using oj L 



II. The second boundary condition is obtained in a similar way with the .result 
being : 


d6(e) 

de 




I L P 
P 


at 


e = 1 


(42) 
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After substituting equation (37) into equation (41) with e = 0 and in (42) 
with e = 1 , the following two simultaneous equations are found: 


• I 

p 


(43) 


*3- p L Ip ) *1 

/ I 63 cos 8 \ 

»3^‘" «3- '1 " “ 


(44) 


Obtaining Nontrivial Solutions 


The procedure for finding non-trivial solutions for Aj and Bj is 
identical to that described for the lateral case* It is straightforward to write 
the expression for the determinant of the 2x2 coefficient matrix in 

2 

equations (43)and (44)* After removing the common 3^ terra (which corresponds to 

the rigid- body mode) from the expression for the determinant, and equating the 
expression to zero, the resulting transcendental equation can be solved by using 
any nonlinear root-finding algorithm* 


CONCLUDING REMARKS 

Equations of motion were derived for free vibrations of a uniform free-free 
beam with a rotational mass attached to each extremity* The simplified equations 
of motion considered in this paper were acquired by assuming the governing beam 
equations to be uncoupled* Using appropriate boundary conditions, nonlinear 
algebraic equations were obtained, the solutions of which yield the modal frequen- 
cies and the mode shapes for lateral as well as torsional vibrations. Computer 
programs were developed for computing the modal frequencies, mode shapes and node 
points* These programs can be used to generate modal data for any number of 
modes. While performing control systems studies based on structural models 
generated, using the computer program, care must be taken to normalize the mode- 
shapes by dividing them by the appropriate L 2 ”norm (i*e*, square-root of the 
Integral of mode-shape squared, integrated from 0 to L over the space 
variable). Mode slopes, which are necessary for control systems studies with 
torque actuators and attitude sensors, can be obtained simply by differentiating 
the corresponding normalized mode shapes with respect to the space variable* 
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WRITE(5»10) FREQ 

— TlI~F0RMAT (18HNATURAL FR ECUENC Y-/ f 1 3. Bi 1X2HH? » / / ) 

C N«AY»N»NRHS AND IFAC ARE ALL VARIABLES USED IN GELIH 

RMAX- 
N-3 
KPHS- 


lPAC-0 


■'a is THE'NORMilllATinN COEFF ICIENT N EEDED TO SOLVE 

FOR THE HOPE SH APES 

A«l. 

S( »_) AND T(_»_) APE THE~~NU HBERS ~IN'THE~ COT FF rC^ 

matrix to be solved' in- GE LIK, 

S(l»l )— Ml*BTA(l)/flU 

S(l/2)-l. 

S(l»3 )— m1*PTA(1)7H0" 

s(2a)— 1. ' 

S(2»2)-I1*BTA{1)**3/(KU*L**2J 

S(2»3)-l. — - - - 

S(3»l )•( 12*9TA (1»**3/ (mj*LA*2)*SIN (BTA'UTFCOS (BTA^DT) 
S(3»2)-(-I2*BTA(l)**3/(HUA[**2T*CO$H (BTA (13)*SINH(B TAmT r 
s (3»3 ) -(-i 2 *PTA (1 )♦ ‘37 {^u*r**2T»5iT5FrnrnrrrT7+cojflmTnnT 

T(lfl)-A ■ ■_■ ■■ _ _ ■ _■ 

T(2»D— I1*BTA(1)**3/(MU*L**2)*A 

T( 3»l )»(I2*BTA(1)**3/(HU*L**2)*C0$ (BTAMD+SINCBTAJl) 

CALL GELlH(NMAy»N,S»NRHS;T;iPIVOT»IFAC#WK,IERR) 


■ BK AND D» ALOHG'VITH'ATRE' THE' CDFFTTIEfTTS^F tMt 
HOPE SHAPE EQUATION 

B-T(l,l) ■_ _ _ _ _ 

c-T( 2 »i ) ■ ■ 

P-T(3»l) ^ ^ IlLZ" 

WRITE(5»15) A»B»C»P»BTA(1) 

FORMATOAHTHE equation of THE MOPE SHAPE ISf//lHP(X/ 
1IN(BETA*X/L)*F15.7»15H*C0S tBETA"*X/L)+/F15.7»16H*SINH 
2.7tl5H*C0SH(BETA*X/L)/7llHWHERE BETA-F1A.7// ) ■ 

JJ-1 _ 

PLTX IS THE increment ADOEP TO'EPSLN'TO' OBTAIN 
■POINTS FOR THE PLOT OF THE HOPE SHAPE 
PLTX-.OOl 


(BETAf-X/D + Fir 


IEC»N,KX»KY»XHIN»XMAX»YNIN»YHAXiPCTPTS»NXHC»NynCi_ 
■ ISYM»Sx;SY,Xj^FF ANP YOFF ARE 'ALL' VARIABLES 

■■' IN INFOPLT. _ 1 

-■'IEC-1 ■ - . -- - 

— N-lOOl 

— - 

— KY-1 I”"”’ ' - 

XHIN-O.O ■ 

XHAX-l. • ■ _ ■ 

- YHIN.-.l 

■'YMAX-.l ~ 

PCTPTS-O.OO ■ ■ - - - 

NXMC-S 


> 

I 


“nymc«V 

ISYH-O — 

'SX-7. 

SY-5. ■' 

■XOFF-,75 — 

Y0FF-.75 

' ~ FPSLN(JJ)~IS~THE~HDRI7o>>'TAL CORRH^ E NT or~THE PLOT 

' WHERE 45 P (K) ' IS~7Hr VEPTKA L . 

E P SLM < 1 ) ■ 0. 

MK)°a*S 1 N( 0 TA( 1 )* E PS L R ( J J) T* B *rt)STST mix? PS CR (T3T)Tr^TrHHWTrg' 

1*EPSLN( JJ) )*0*COSH(B TA (11*EPSCN'( JJTr~ 

'IF (K.EO.l) ‘GO T0_5 

the following I F'~S mEHENT^S~Tj~S' E~0 ~rcrQ'6TnN A 

FIRST GUESS' AS TO ThE LOCATION 0'P~~tHE NO DES AND 

— SURPOUTINE NOOE 'USES' THE NEBTOIT ITAP'H5aH~P gTFDD TO 

— OBTAIN THF'TRUE'POOTS; 

“ IF (R ( K ) ♦ B ( K-1 ) . LE^OnTA rr-NODETrPSLKTJTriTAl 

“EPSLN { JJ + 1 )»EPSLRr j'Jl»bLTX 

~JJ-JJ+1 ■ ■ 

"CONTINUE 

'CALL I N F 0 P L T ( I E C » M » FPSOnTnnrrirriTrKY7yRITT7^yTX7TyT KlYHAX» 

TPCTPTS»NXHC» ?HEPSLN»KY"C7 AHR TF)» rS7M>SX7 SY7XOFFrYOFF) 

“RETURN 


V 


SUBP OUTI NE“FVACf 0TA7 XYT1 [ 

"CO** MON/ VALUES/ I1/IZH/HI»R?7T'U7FI» ATBTCVOT^D 

-REAL T1»IZ»L>H1»R2»HU ' ' 

dimension BTAaifXYZm 

-V.fP»rn*;tftTAm)*COSH(BTAnn-21*L**A *RU*PA>-<I2»BTA-m*CDS^( BT- 

“I 4 { 1 ) ) *SINH(BTAmi-Z*PTA(l 1‘COSFrBTAm )*SlNtBi Ati» » I'LATAPr z*( 

~Z~ 2*BTAm*C0S(BT4(n )*STKH(BTA(l))-2*BTA(n*C0SHIBTAm)*SIN(BTA 

"3 — (1 ) ) ) ♦L**A*M1 + 1 12* (-2*BTA(rV*P3*COSCBTA (ITV*SINH (BTA (1) )-2*BTAT 

“A — l)**3*C0SH(BTA{l))*SIK(PrA(I)))*Il*C-2ABTAm**3*C0S(BTA(in*^ 
"5 — NH(PTA(in-2*BTA(n**3*COSH(BTA(l))*SIN(BTA(in nn**2)*MU**3 

— Y»{ (( I2*(-2*BTA(1)**A*C0S(6TAM) ) pCOSHC BTA ( I ) ) -2*BT A ( 1 ) **A ) - AXI^X^ 
"7 T A ( 1 ) ♦ ♦ A * c 0 ST 6 Trm'l^CTISnTFTrrTTTTTPf *Z=TTFTAU)* ♦ 2 ♦SlNTBTmTT 

-ft — ♦siNH( 6 TA(l))*L**A*MlJ*M 2 +(ll*(“ 2 *BTA(l)**A*COS(BTA(in*COSH(BT 
-S ■ A(in‘2*5TA(l)**A)-Ail2*BTA(l)**A*COS(PTA(in*COSH(BTA(in)*L** 

I 2***l + A*Il*I2*BTA(l)**6*SlN(BTAm)*SINHtRTA(im*MU**2 

-;^’*i!TMl)P*5*COSH(BTAm)*SlN(BTA(l))-2*PTA(l)P*5*COS(BTA(in*SIN 
H( BTA(1) J )*Il*(2*PTArrV»*5*C0SHTBTA m ) ASIfUPTA (m-2*BTAm 
*COSIBTAm)*SINH(BTAm) ) )*L*P 2 *M 1 + I1*I2* «2*BTA(1)**7*C0S(BTA1 

V 1) )*SINH(8TA(1))*2*PTA(1)p*7*C0SH(BTA(1))pSIN(PTA(1) ) ))*M2+I1*I 



2* (2*0T4(1 )**7*C0S(BTA(1) ) ♦SINH ( B TA ( 1 ) ) *2*61 A ( 1 ) ♦♦7 *CDSH ( BT A ( I ) 

— )*siN(BTA(i j ) )*f<r)*»u+n*i2*T7^Errrrn>*'B»r'03‘riTrn:)T^C7STTrmTi 

))-2*BTA(l)**9)*Hl*H2 

XY7 (1 » ■ ( X ♦y+ 7 r/ ( L**A*HU*A^J 
PETU9N 
’ENO'^' 


_ SUB<li 

_corno 

REAL 
DIME' 
- U 
T“ L* 
~ 2 ~~ 2 * 
~3 N( 

1 ) 

-fc - )) 
“V-' ( 

"7 — ♦* 

-f, 1 , 

"3 2« 

■|— (1 
2 - 

-< — 

'• — CG 

“> A( 

~7 )♦ 

-“V“- ( 

S SI 

T" a 
“2 — *i 
-3 BT 
-Jf 

-5 )♦ 

~6" H( 
- 7 — ♦! 
"9 ' S( 
~ DPPIM 
■“RtTUR 
“END 


OUTINE JAC0«(BTA»DPRJME1 

N/VALUES/ Il>I 2VL;Hi;h2^FU7n >A»B»C/D>CPD ~ 

Il/I2»r;>U»>''2^nU 

NSION BTA(l)iDPRiP£(l7T) ^ 

•( 2 *COS(BTA(in*SINMl BT Anr FI»‘CT5 S HT BT ATTITTS Tm BTXTlTrT * 
*4*f*(j**^V( ( (2*CnS (0 TaT1T)-A*BTa (1 JTEsfRTBf Afl JTTASTTTHfBT'An ) )- 
cosh(8ta(i) )*siN(BTAmTrnr^*^H2vn2*cosTBTi(r))-A>6TA(rr<^r 
PTA(in)*SINH(BTA(nF2*C0St7rRTA(in*SIfjrBTA'(llTr*L'«A + Nl+T(- 
1 2=6 * 1 1 ) ♦ P T i'( IT^tstctJBTB T A'( IT) ’TT R h (BTAU)) + (- fc*l 2 -f' *TIT* B T A T 
**2*C0SH(BTAm)*SIN(BTA'nT)+T-«*T2=A*rn*BTAnT**3»C0SlBTAri 
*COSH(BTAm) J*L**2J*hU**3 
( ( (“8*PTA(l)*SIN(eTA (in=A*BTAri) 

2*COS(8TA(l) ) )*SIHH(BTA(r) )-A*RTArn**2*CDSK(BTATm*SIN(6TAr 
))»L**A*N14( (-2*1 2=A*ir)*PTA (n **'4*C0S (PTA (in*SINH(BTA I m+r 
1 2 ♦ A * 1 1 j ♦ p T A »c n smTnrrrrrriTTTTTriTT 4T=F» r2-i £ *tittbt7t 

)**3*C0S(BTA(1) )*COSH(BTA(l))=B*T2*BTAm**3)*L* + 2)*M2 + ( (-A*r 
2*IU*aTA(l)**A*COS(BTA(l))ASINH(0TA(l)) + (A*T2'+2*Il ) + BTA(l 
CnSH(BTA(l))*SIN'(BT A (!)) + (-! 6*1 2-e*Il)*BTAm**3*COS(BTAm7*“ 
SH(BTA(1) )-P*Il*BTA(l)**3)*L**2*Nl+ (2A*T1*T2+BTA(1 )**5^SINIBr 
l))4A*Il*I2*BTAri)*«6*CnS(BTAmn*SINH{nTAmH-AAIl*I2*BTArr 
♦ 6 * C D SH X B T A mi ♦ STNXBTlTmT'TD » * 2 
(({(A*I24A*Il)*BTAm*«5*“" ' ■ 

N(RTA(1))*(-10*I2-10*I1)*BTA(1)»*A+C0S(BTA(1) n*SINH(BTA{l))+ 
0 + I2 + 10*n)*BTA(l)**A*C0SH(PTA(l))*SINtBTA(l)))*L* + 2*Nl + lA + Il 
2*BTA (1 )**f*COS(BTA (1) )*SINH(BTA (1 ) )+lA*Il*I2*BTA(l)**6*COSH( 
A(1) )*SIN(BTA(1) )4A*Il*I2*0TA(l)*A7*COS(RTA(in*COSH(BTA(l )) T 

2 4 » 1 A * 1 1 * 1 2 * BTa (1 ) * * 6* C DS f BTATmTf’rrNPTB TTnTTnTA^' 1 1 ♦ 1 2 TTEF 

»6*C0SH(BTA(1))*SIN(PTA(1 n4A*Il*I2*RTA<n**7*C0S{6TA(l))*C0S 
BTA ( 1 ) ) ) *N1 ) *“U4 (2*I1*I2*BTA(1)*A0*COS(BTAC1))*S1NH(BTA(1))-2 
1*I2*6TA(J )**fl*COSM{BTA(l))*SIN(BTAn ) ) 4l6* II* I2*B^(1 ) ♦*7*C0 

BTA XI ) )*C0SH{BTA(H )-16*Il*I2*BTA (1)**7J*N1*B2 — — 

E(l4l)«(U*V4W)/(L**A*MU**A)' 


— SUBROUTINE NODE X EPSLN» JJ» BTA ) “ ' 

— C0MH0N/VAlUES/Il»I2»L»Nl»N2»riU»El> A>B>C»D4CR0 


ABLE ANO'BAKEP APE ' THF EXTERNALS" THAT C ONTATN-THE' 
NODE SHAPE EQUATION AND ITS DERIVATIVE WHICH ARE" 



> 

I 

ON 


c Used by nonlin# “ 

' EXTERNAL ABLE#BAKER 

OI^^ENSION BTA (l)#BPSLN(2bOO)#ANDb( l)»WK (4)#AA(1#1)#CRD(1) 

ZLpEAL ^ 

ANODm-EPSLN(JJ) 


“C NMAX,N,NSIG/IKJAC AN D“n N J5T S E CT‘6 Y~NO N Cl N"i 

NmAX «1 * - ' . 

N-1 

' NSTG-14 ■ 

lNjAC-1 

, ITMAX -10000 

^ ■ CRD (1 )-BTA(l) _ 

"CALL NONLIN(NHAX#N,ANOD/N$lG#iNJ'AC,BAKER/ITK'AX^'ABLE^Wkf AA'^IERn' 
"X-ANnP(l)*L 
WRITE(5#10) X 

I'O FQ0MAT(23HA NODE' I'S L0C4TE0“AT“X»F l4T5 J 
■ ' RETURN 


‘ SUBROUTINE AalFl anooVrr") 

CCHhON/VALUES/ 1 1 » 1 2 i*L7 M IJV. Z7MU7E ITA7ET# Z,d,ZPV 
DIMENSION ANOO(l)^RR(l )fCRDCn 

RR (l)-A*SIN(CRD(l)*ANOOm )+B*CDSrCRD(r)>ANDD(rTT^C^3 TNH(i:RDTr) *' 

1AN0D( 1) )+n*COSH(CRD(l)*ANODm ) 

RETURN 

END 


■ SUBROUTINE PAKFR ( ANOD# ORR J “ 

~ C 0 HH 0 N/VALUES/Il#I 2 #L>Nl>M 27 MU 7 HrrA 7 B 7 CT 0 #CRD 
■■■ DIMENSION ANOO(l)#DRP(lVinCPDfl) 

■ DRR(Ul)-A*CRD(l)>CnS(CRDTTT»'WnD(l) )-B*CROll)*SIN(CRDriT*’A-NDtrm 

l)+C*CPD(l)^COSH(CRDm*ANODm)+DKRD(l}*SINH(CRO{l)+ANOD(l) ) 

' RETURN "■ “ 

' ENO ^ 




1 1 


r 


iEXAMPLE^P‘‘fa \TERA l7~V^fil^I ^ PROGR AM- 


C _ THE FOLLOWING 0 VARIABLES ARE THE INPUT VARIABLES. 

^19 C _ ^ ---.-. 

C HU- HASS OF THE BEAM ALONE 

' L-LENGTH OF THE BEAN 

C „ Nl-NASS AT X-0. END OF BEAM 

_ C N2-NASS AT X-L END OF BEAN 

JO C . n-INERTIA OF Ml 

C I2-IHERTIA OF M2 ’ 

C _ EI-STIFFNESS OF THE BEAM _ 

_ C NMOOE-THE NUMBER OF MOOES SOLUTIONS ARE NEEDED FOR» STARTING 

C _ WITH II AND ENDING WITH INMOOE. 



MU-177. " 

L-100. “ 

Mi-eoooo.o 

M2-250. "" 

_30 Il-93flOU1.7 _ *_ ; ' 

. 12-250. ; ' ' ' 

El-20000000.0 ' 

NMr.DE-5 ‘‘ “ " ^ 


> 

'U 

n) 

13 

G- 

H* 

X 


tsi 

• 


w 

CD 


3 


TJ 




T3 

C 


THES 

_BEAM 

THE 

_THE 

_the 

-.the 

_TWE 

_THE 

THE 


f JRB_THE_50LUTI0NS_FCR LATERAL VI BR ATIOM^OL JHE 
.WITH THE FOLLOWING CHARACTERISTICS. 


MASS OF THE BEAM - 
LENGTH OF THE BEAM- 
MASS LOCATED AT X-0. - 
MASS LOCATED AT X-L - 
INERTIA_0F_THE MASS^AT_ 
inertia of the mass AT X-L 


177.000 

100.000 

eoooo.ooo 

. 250.000 

X-0. - _ 93801Al.700_ 
250.000 


STIFFNESS OF THE BEAM • 20000000.000 


(D 

G 

G. 

O 

c 

n 

c 


o 

►TJ 

w 

m 


^ - 


_R1S U L T S _F OR M CO E _ • __1 


i j<atur'al_Meo'uen cY- .oooooboo hz 


i JTHE JQUATION OF THE MO'DE S*PE ISi ^1' 

f 

l~R ( X/ L)» ~~ 1,00 00000»S1N( BETA»X/l)» ^.000600i»CdS~( 8ETA*x7n» 






J HH. J OU AT I OH ;'OF^ Til E j) 0? £ J” * ^ . I S • 

gU/L)» I.OOOOOOO-^SIHCBETA^X/ D^ ^r.34‘?229S»Cbr('BETl '^ 

7l.2622367‘>SINH(SETA*X/L)y"_l 1.3 4521 A5*C OSH (B ETA^l/Li 

7hE re be T A- 1 2 46 94 2 6 7 ' y 


a'NOO'E TS ibCATEb' at "X-" ■■i7.'79221 


,R E s ii L t s, . F OR . noo E Jl .11 


NATURAL Ff TEO UENCY- .e 3'667 199'HZ~ 


.the! EbbAtjONloF . the_lvo60hape7is I 

"r (x7 Li -“I iTbo6bdb6*slN ( beta ♦x7lT»~I ^6 6 2 b 'm ♦ c os ( b e t a * x / l m 

»1.005 B983*SIHH( B ETA *X/L )» i.6 009016*C OSH ,( BET A*X/ L> 

!V her V BE TA* 4 .Idoi 594 1 


A_NbOE_IS LOCATEO AT X- I .9i457“ 

A NODE IS LOCATED At ~X* 96 .441 20 


RESULT S. FOR ..MOd E_L..A 


natural' FREQUENC Y* ~ 2.66645 526 HZ 


L THE EQUATION OF THE~MODE SHAPE ISt 


B-2 







A^NOOE^IS LOCATED ATJ(- •361A0 

A NODE IS^LQCATEO^AT X« 53.64698 

A'mOOE is located Af X» 90.77667 



MAfURAL FREQUENCY- 5.52140730 HY 


J H A TTo ?L Of ^Th_e . m od e . s h a p eJ Y s'j 

117/ U -3 '“l . OOOOOOO^SINCBETA^X /C V^ «>1.00T7B34^Y(BTf r»^^ 

ri*-P0Q36p.l#SLNH^ ^l-p003>77^C0SH(iElAlXn 

VhER E J rTA^l“ 1671590663 


A NODE IS LOCATED AT X- 
A_NOOE IS LOCATEO_AT X- 

A node'is"lqcated at X- 

A NODE IS_L0CATED1atJ^X- 



.20791 

30.79191" 

69.62796 

99.45123 












G m 



Figure 7.- Output from FFBEAM plotting displacement R(E) versus 
nondimensional position variable E for third mode. 
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APPENDIX C. 


MACSYMA PROCEDURE 

This appendix lists and explains the code used to derive the equations 
necessary to solve for the eigenvalues for free lateral vibration of the configura- 
tions described earlier. The following is a list of a batch file which can be run 
on MACSYMA to duplicate the derivation that was used for lateral vibrations: 

1. A*SIN(BETA*E)+B*COS(BETA*E)+C*SINH(BETA*E)+D*COSH(BETA*E); 

2. DIFF(D2,E,1); 

3. DIFF(D2,E,2); 

4. DIFF(D2,E,3); 

5. D5-M1/MU*BETA**4*D2=0.; 

6. SUBST(0. ,E,D6); 

7. D5+M2/MU*BETA**4*D2=0.; 

8. SUBSTd . ,E,D8) 

9. D4+I1/(MU*L**2)*BETA**4*D3=0.; 

10. SUBST(0. ,E,D10); 

11. D4-I2/(MU*L**2)*BETA**4*D3=0.; 

12. SUBSTd . ,E,D12); 

13. C0EFMATRIX([D7,D9,D11,D13] ,[A,B,C,D]); 

14. DETERMINANT^ %); 

15. RATSIMP(%); 

16. DET=D16; 

17. DIFF(D16,BETA,1); 

18. RATSIMP(%); 

19. DPRIME=D19; 

20. F0RTRAN(D17) ; 

21. FORTRAN(D20); 


The following is an explanation of each line individually: 
Line Number Definition 


1 

This 

line 

2 

This 

line 

3 

This 

line 

4 

This 

line 

5 

This 

line 

6 

This 

line 

7 

This 

line 

8 

This 

line 

9 

This 

line 

10 

This 

line 

11 

This 

line 

12 

This 

line 

13 

This 

line 

14 

This 

line 

15 

This 

line 

16 

This 

line 


simply writes equation (8). 

differentiates equation (8) with respect to e. 
differentiates equation (8) twice with respect to e. 
differentiates equation (8) three times with respect to e. 
writes the first boundary condition: equation (16). 
substitutes zero for e to satisfy the boundary conditions, 
writes the second boundary condition: equation (18). 
substitutes one for e to satisfy second boundary condition, 
writes the third boundary condition: equation (22). 

substitutes zero for e to satisfy third boundary condition, 
writes the fourth boundary condition: equation (23). 

substitutes one for e to satisfy fourth boundary condition, 
generates the four by four matrix, 
generates the determinant of the matrix, 
simplifies the determinant, 
writes the determinant in the form DET=.... 
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Line Number 


Definition 


17 This line differentiates the determinant with respect to 3. This is 

needed to use the Newton-Raphson method of finding roots. 

18 This line simplifies the derivative. 

19 This line writes the derivative in the form DPRIME=.... 

20 This line writes the determinant in FORTRAN code. 

21 This lines writes the derivative in FORTRAN code. 


This file can be run on MACSYMA by typing the command; 

BATCH( DKROBE , LATRAL , DSK , LRC) ; 

This must be done as the first command in a session: on line ”C1'*. 
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Figure !•- A differential beam element showing shears and moments resulting 
from displacement V(x,t). 


Figure 2.- (a) The beam configuration considered in this paper with 

rotating masses m^ and m 2 at either end. Parts (b) and 
(c) show the force and moment reactions at the ends. 
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